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Abstract 



We confront a very wide body of HERA diffractive electroproduction data 
with the predictions of the colour dipole model. In doing so we focus upon 
three different parameterisations of the dipole scattering cross-section, whose 
parameters are fixed by analysis of DIS structure function data. This analysis 
strongly suggests the presence of saturation effects, but is not definitive because 
the conclusion requires the inclusion of data at low-Q 2 values. 

Having fixed the parameters of the models from the DIS structure function 
data, the resulting dipole cross sections can be used to make genuine predic- 
tions for other reactions. Good agreement is obtained for all observables, as is 
illustrated here for deeply virtual Compton scattering(DVCS) and diffractive 
deep inelastic scattering(DDIS). There can be no doubting the success of the 
dipole scattering approach and more precise observations are needed in order 
to expose its limitations 



1 Introduction 

In the colour dipole model [cH 0] , the forward amplitude for virtual Compton scatter- 
ing is assumed to be dominated by the mechanism illustrated in FigO]in which the 
photon fluctuates into a qq pair of fixed transverse separation r and the quark carries 
a fraction z of the incoming photon light-cone energy. Using the Optical Theorem, 
this leads to 



*This talk is based on work carried out in collaboration with J.R. Forshaw and 
R. Sandapen, [HE] 





for the total virtual photon-proton cross-section, where \I/^' T are the appropriate 
spin-averaged light-cone wavefunctions of the photon and a(s*, r) is the dipole cross- 
section. The dipole cross-section is usually assumed to be independent of z, and is 
parameterised in terms of an energy variable s* which depends upon the model. 




Figure 1: The colour dipole model for elastic Compton scattering "f*p — > 7*p. 

Thus using Eq.((l]) we are able to compute the deep inelastic structure function 
F 2 (x, Q 2 ). The power of the dipole model formulation lies in the fact that the same 
dipole cross-section a(s*, r) appears in a variety of other observables which involve 
the scattering of a real or virtual photon off a hadronic (or nuclear) target at high 
centre-of-mass (CM) energy. The largeness of the CM energy guarantees the factor- 
ization of scattering amplitudes into a product of wavefunctions and a universal dipole 
cross-section. In this paper we wish to test the universality of the dipole cross-section 
using a wide range of high quality data collected at the HERA ep collider. Moreover, 
we also wish to examine the extent to which to data are able to inform us of the role, 
if any, played by non-linear saturation dynamics. 



2 The dipole cross-section 

We now turn to the three different models used to describe the dipole cross-section. 
Before doing so however, we shall first discuss our choice of photon wavefunction. 
For small r, the light-cone photon wavefunctions are given by the tree level QED 
expressions 

These wavefunctions decay exponentially at large r, with typical r-values of order Q^ 1 
at large Q 2 and of order mj 1 at Q 2 = 0. However for large dipoles r 1 fm, which are 
important at low Q 2 , a perturbative treatment is not really appropriate. In this region 
some authors j3] modify the perturbative wavefunction by an enhancement factor 
motivated by generalised vector dominance (GVD) ideas [HI E], while others P 
achieve a similar but broader enhancement by varying the quark mass^. In practice pQ, 
the difference between these two approaches only becomes important when analysing 

TFor a review of the colour dipole model, including a fuller discussion of these points see Forshaw 
and Shaw (TOj - 
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the precise real photoabsorption data from fixed-target experiments jTTj. Since we will 
not consider these data here, we will adopt the simpler practise of using a perturbative 
wavefunction at all r- values, and adjusting the quark mass to fit the data. 

Turning now to the dipole cross-section, all three models are consistent with the 
physics of colour transparency for small dipoles and exhibit soft hadronic behaviour 
for large dipoles. As stated above, the model parameters are determined by fitting 
only to the DIS structure function data. Since the details of all three models have 
been published elsewhere, we shall here summarise their properties only rather briefly. 

2.1 The FS04 Regge model |TJ 

This simple model, due to Forshaw and Shaw jT], combines colour transparency for 
small dipoles r < r with "soft pomeron" behaviour for large dipoles r > r\ by 
assuming 

a(x m ,r) = Anr 2 x^ H for r < r and 

= A s x^ Xs for r > r 1} (2) 

where x m = (Q 2 + 4m 2 ) /(Q 2 + W 2 ) . For light quark dipoles, the quark mass rrif is a 
parameter in the fit, whilst for charm quark dipoles the mass is fixed at 1.4 GeV. In 
the intermediate region ro < r < r±, the dipole cross-section is given by interpolating 
linearly between the two forms of Eq. (J2J) . 

If the boundary parameters r and r\ are kept constant then this parameterisation 
reduces to a sum of two powers, as might be predicted in a two pomeron approach, and 
can be thought of as an update of the original FKS Regge model j^] to accommodate 
the latest data. It is plainly unsaturated, in that the dipole cross-section obtained 
at small r- values and fixed Q 2 grows rapidly with increasing s (or equivalently with 
decreasing x) without damping of any kind. 

2.2 The FS04 Saturation model p] 

Saturation can be introduced into the above model by adopting a method previously 
utilised in jT2J. Instead of taking tq to be constant, it is fixed to be the value at which 
the hard component is some specified fraction of the soft component, i.e. 

v{x m , r )/a(x m , n) = / (3) 

and / instead of tq is treated as a parameter in the fit. This introduces no new 
parameters compared to the Regge model. However, the scale r now moves to lower 
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values as x decreases, and the rapid growth of the dipole cross-section at a fixed, 
small value of r begins to be damped as soon as r becomes smaller than r. In this 
sense we model saturation, albeit crudely, with r the saturation radius. 

2.3 The CGC saturation model [TH] 

In addition we shall consider the CGC dipole model originally presented by Iancu, 
Itakura and Munier [T^j. This model aims to include the main features of the "Colour 
Glass Condensate" regime, and can be thought of as a more sophisticated version of 
the original "Saturation Model" of Golec-Biernat and Wusthoff [5]. The original 
Iancu et al dipole cross-section was obtained using a three flavour fit to the DIS data 
Here we use a new four-flavour CGC fit due to Kowalski, Motyka and Watt |14j . 



3 Structure function data 

The parameters of the FS04 models were determined by fitting the recent ZEUS F 2 
data ^B] in the kinematic range 

0.045 GeV 2 < Q 2 < 45 GeV 2 x < 0.01 (4) 

whilst the CGC fit of [Hj is to data with Q 2 > 0.25 GeV 2 (the other limits are as 
for FS04). The corresponding HI data [IB] could also be used, but it would then be 
necessary to float the relative normalisation of the two data sets. We do not do this 
since the ZEUS data alone suffice. The resulting parameter values are tabulated in 
the original papers; we do not reproduce them here, but confine ourselves to some 
general comments. 

The best fit obtained with the FS2004 Regge model is shown as the dashed line in 
Figure El (left). As can be seen, the quality of the fit is not good, corresponding to a 
X 2 /data point of 428/156. This is not just a failing of this particular parameterisation. 
We have attempted to fit the data with many other Regge inspired models, including 
the original FKS parameterizations 011^], without success. 

A possible reason for this failing is suggested by the solid curve in Figure El (left), 
which shows a the result of fitting only to data in the x-range 5 x 10 -4 < x < 10 -2 , 
and then extrapolating the fit to lower x values, corresponding to higher energies at 
fixed Q 2 . As can be seen, this leads to a much steeper dependence at these lower 
x-values(i.e. at higher energies) than is allowed by the data at all Q 2 . An obvious way 
to solve this problem is by introducing saturation at high energies, to dampen this 
rise. This is confirmed by fitting with the the FS04 saturation model, which yields a 
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Figure 2: Comparison of our new dipole model fits to a subset of DIS data [T^j. Left: 
No saturation fits. FS2004 Regge dipole fit (dashed line) and (solid line) a fit of the 
same model to data in the restricted range 5 x 10 -4 < x < 10~ 2 , extrapolated over 
the whole x-range x < 0.01. Right: Saturation fits. FS2004 saturation fit (solid line) 
and the CGC dipole model (dot-dashed line) 

X 2 /data point of 155/156. The resulting fit is shown in Figure El (right), which also 
shows the very similar results previously obtained using the more sophisticated CGC 
model. 

It is clear from these results that the introduction of saturation into the model im- 
mediately removes the tension between the soft and hard components which is so 
disfavoured by the data. However, it is important to note that this conclusion relies 
on the inclusion of the data in the low Q 2 region: both the Regge and saturation 
models yield satisfactory fits if we restrict to Q 2 > 2 GeV 2 , with % 2 /data point values 
of 78/86 and 68/86 respectively. 

At this point we have three well-determined parameterisations of the colour dipole 
cross-section, which have been used by Forshaw, Sandapen and Shaw [2] to yield 
predictions for other processes. In the next sections we shall take a look at the results 
for Deeply Virtual Compton Scattering (DVCS) and the diffractive structure function 
i^P^t .We always choose to show the Regge fit, even though it does not fit the F 2 data 
particularly well, in order to indicate the discriminatory power of the data. We stress 
that in all cases, the photon wavefunctions and dipole cross-sections are precisely 
those determined from the fits to F% data, without any adjustment of parameters. 
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Figure 3: Comparison of the ZEUS DVCS data [21] with the predictions of the three 
models discussed in the text. Left : Q 2 dependence at W = 89 GeV. Right: W 
dependence at Q 2 = 9.6 GeV 2 . 



4 Deeply virtual Compton scattering 

In deeply virtual Compton scattering, 7*+p — > 7+p, the final state particle is a real as 
opposed to a virtual photon and dipole models provide predictions for the imaginary 
part of the forward amplitude with no adjustable parameters beyond those used to 
describe DIS. To calculate the forward cross-section a correction for the contribution 
of the real part of the amplitude has to be included. This correction was estimated in 
[TH] and found to be less than « 10% and of a similar size in different dipole models. 
Predictions for the measured total cross-sections are then obtained using 

(5) 

t=o 

where the value of the slope parameter B is taken from experiment^. 

The predictions of all three dipole models are compared with the ZEUS data j2I] in 
Fig0", where we take a fixed value B = 4 GeV -2 which is compatible with their 
data. Bearing in mind this normalisation uncertainty, the agreement is good for all 
three models, although significant differences between the models appear when the 
predictions are extrapolated to high enough energies, as one would expect. Similarly 
good agreement is found for the HI data (22j (see [2]). 

* Predictions for vector meson production, which require a discussion of the vector meson wave- 
functions, can be found in [181 12] . 

§For an alternative investigation of the link between low- a; DIS and DVCS and other exclusive 
proceses at high energy, see Kuroda and Schildknccht [20] , 

^Note that throughout this paper the curves labelled 'FS04 no sat' correspond to the predictions 
of the FS04 Regge model. 



<7l,t(j*P ~> TP) 



1_ da T < L 
B dt 
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5 Diffractive deep inelastic scattering (DDIS) 



WW 




Figure 4: The qq and qqg contributions to DDIS. 

To conclude our study we turn to the diffractive deep inelastic scattering (DDIS) 
process 

where the hadronic state X is separated from the proton by a rapidity gap. In this 
process, in addition to the usual variables x and Q 2 there is a third variable M x . In 
practice, x and M x are often replaced by the variables xjp and j3: 

M 2 X + Q 2 jc_ Q 2 

X]P ^ W 2 + Q2 P- Xp ^ M 2 X + Q 2 • W 

In the diffractive limit s > Q 2 , m 2 x and so xp < 1. 

In the dipole model, the contribution due to quark- ant iquark dipoles to the structure 
function can be obtained from a momentum space treatment as described in 

|231 IH]. However, if we are to confront the data at low values of j3, corresponding 
to large invariant masses M x , it is necessary also to include a contribution from the 
higher Fock state qqg. We can estimate this contribution using an effective "two-gluon 
dipole" approximation due to Wiisthoff [221, as illustrated in FigJU 

Again, the predictions obtained in this way involve no adjustment of the dipole cross- 
sections and photon wavefunctions used to describe the F2 data. We are however free 
to adjust the forward slope for inclusive diffraction (B) within the range acceptable 
to experiment, which means that the overall normalisation, but not the energy de- 
pendence, of is free to vary somewhat. We take B = 6.8 GeV~ 2 when making 
our CGC and FS04 saturation predictions and B = 8.0 GeV -2 when making the 
FS04 Regge predictions. Note that a value of 8.0 GeV -2 is rather high compared to 
the « 6 GeV -2 favoured by the HI FPS data [21] although it is in the range allowed 
by the ZEUS LPS data [213 • The need for a larger value of B for the FS04 Regge 
model arises since the corresponding dipole cross-section is significantly larger than 
the FS04 saturation model at large values of r and this enhancment is magnified in 
inclusive diffraction since it is sensitive to the square of the dipole cross-section. We 
should also bear in mind that the tagged proton data are subject to an overall ~ 10% 
normalisation uncertainty. We are also somewhat free to vary the value of a s used to 
define the normalisation of the model dependent qqg component, which is important 
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Figure 5: Model predictions compared to the ZEUS LPS data [23] . 

at low values of (3. Rather arbitrarily we take a s — 0.1 and take the view that the 
theory curves are less certain in the low (3 region. 

In FigJSJwe compare the recent ZEUS LPS data [53] on the xp dependence of the 
structure function F^^ 1 at various fixed Q 2 and (3 with the models" . The agreement 
is good except at the larger x p values. Indeed, the x 2 values per data point are very 
close to unity for all three models for xp < 0.01. Disagreement at larger xp is to 
be expected since this is the region where we anticipate a significant non-diffractive 
contribution which is absent in the dipole model prediction. Note that the three 
models produce similar predictions at larger values of j3. 

Comparison to the HI data with tagged protons [21] is to be found in Fig|U]and Figd 
The story is similar to that for the ZEUS data and the evidence for an overshoot of 
the CGC and Regge model predictions at low (3 is strengthened. The Q 2 = 2.7 GeV 2 
panes in Fig|7| illustrate this point the best. Again, we should not interpret this as 
evidence against these dipole models due to the uncertainty in the qqg contribution 
in the low (3 region. The agreement between all models and the data at larger values 
of (3 and low enough xp is satisfactory.**. 

"Predictions for the original Iancu et al CGC model have previously been published in |26| . 
**A fuller discussion, including comparison with the ZEUS FPC data [22j and the HI My data 
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Figure 6: Model predictions compared to the HI FPS data: xjp dependence |24j . 

In summary, the DDIS data at large enough (3 > 0.4 and small enough xjp < 0.01 
are consistent with the predictions of all three dipole models. However the data 
themselves would have a much greater power to discriminate between models if the 
forward slope parameter were measured to better accuracy. At smaller values of f3, 
the data clearly reveal the presence of higher mass diffractive states which can be 
estimated via the inclusion of a qqg component in the dipole model calculation under 
the assumption that the three-parton system interacts as a single dipole according 
to the universal dipole cross-section. The theoretical calculation at low (3 must be 
improved before the data in the region can be utilised to disentangle the physics of the 
dipole cross-section. Nevertheless, it is re-assuring to observe the broad agreement 
between theory and data in the low /3 region. 

[2Hj is given in 0. 
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Figure 7: Model predictions compared to the HI FPS data: j3 dependence |24j . 



6 Conclusion 



The dipole scattering approach, when applied to dim-active electroproduction pro- 
cesses, clearly works very well indeed. The HERA data now constitute a large body 
of data which is typically accurate to the 10% level or better, and without exception 
the dipole model is able to explain the data in terms of a single universal dipole scat- 
tering cross-section. Perhaps the most important question to ask of the data is the 
extent to which saturation dynamics is present. Although the F 2 data suggest the 
presence of saturation dynamics the remaining data on exclusive processes and 
on F®^ 1 are unable to distinguish between the models we consider here: these data 
are therefore unable to offer additional information on the possible role of saturation. 
We do note that a more accurate determination of the forward slope parameter in 
diffr active photo/electro-production processes would significantly enhance the impact 
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of the data. However, it is hard to avoid the conclusion that only with more precise 
data or with data out to larger values of the centre-of-mass energy will we have the 
chance to make a definitive statement on the role of saturation without the inclusion 
of the low Q 2 -F 2 (x, Q 2 ) data in the analysis. 
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